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ABSTRACT 

The detection of planar structures within the satellite systems of both the Milky Way 
(MW) and Andromeda (M31) has been reported as being in stark contradiction to 
the predictions of the standard cosmological model (ACDM). Given the ambiguity 
in defining a planar configuration, it is unclear how to interpret the low incidence of 
the MW and M31 planes in ACDM. We investigate the prevalence of satellite planes 
around galactic mass haloes identified in high resolution cosmological simulations. We 
find that planar structures are very common, and that ~10% of ACDM haloes have 
even more prominent planes than those present in the Local Group. While ubiqui¬ 
tous, the planes of satellite galaxies show a large diversity in their properties. This 
precludes using one or two systems as small scale probes of cosmology, since a large 
sample of satellite systems is needed to obtain a good measure of the object-to-object 
variation. This very diversity has been misinterpreted as a discrepancy between the 
satellite planes observed in the Local Group and ACDM predictions. In fact, ~10% of 
ACDM galactic haloes have planes of satellites that are as infrequent as the MW and 
M31 planes. The look-elsewhere effect plays an important role in assessing the detec¬ 
tion significance of satellite planes and accounting for it leads to overestimating the 
significance level by a factor of 30 and 100 for the MW and M31 systems, respectively. 
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1 INTRODUCTION 

While the Universe at large may be homogeneous and 
isotropic, on galactic scale the distribution of galaxies is 
highly anisotropic. This is most readily seen in the spa¬ 
tial and kinematical distribution of the Local Group (LG) 
satellites. In the MW, the 11 “classical” satellites define 
a thin plane (Lynden-Bell 1976) and some of the fainter 
satellites, tidal streams and young globular clusters have 
an anisotropic distribution reminiscent of this plane (Metz, 
Kroupa & Jerjen 2009; Pawlowski, Pfiamm-Altenburg & 
Kroupa 2012). Many members of this “disk of satellites” 
have a common rotation direction and it has been claimed 
that the plane is a rotationally stabilized structure (Metz, 
Kroupa & Libeskind 2008; Pawlowski & Kroupa 2013). Sim¬ 
ilarly, the spatial distribution of satellites around M31 is 
anisotropic (Koch & Grebel 2006; McConnachie & Irwin 
2006), with 15 out of 27 satellites observed by the Pan- 
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Andromeda Archaeological Survey(PAndAS; McConnachie 
et al. 2009) located in a thin plane. Out of the 15 members 
of the plane, 13 of them share the same sense of rotation 
(Ibata et al. 2013, hereafter Ibatal3). 

Anisotropies in the distribution of satellites are a clear 
prediction of the A cold dark matter (ACDM) paradigm 
(Libeskind et al. 2005; Zentner et al. 2005; Libeskind et al. 
2009, 2011; Deason et al. 2011; Wang, Frenk & Cooper 2013; 
Sawala et al. 2014). Such flattened satellite distributions, 
dubbed “great pancakes”, can arise from the infall of satel¬ 
lites along the spine of filaments (Libeskind et al. 2005; 
Buck, Macci’o & Dutton 2015), which in turn determine the 
preferential points at which satellites enter the virial radius 
of the host halo (Libeskind et al. 2011, 2014). The imprint 
of anisotropic accretion is retained in the dynamics of satel¬ 
lites, with a significant population co-rotating with the spin 
of the host halo (Libeskind et al. 2009; Lovell et al. 2011; 
Cautun et al. 2015). 

Although flattened satellite distributions are common 
in ACDM, configurations similar to those of the MW and 
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M31 are infrequent. Wang, Frenk & Cooper (2013) fonnd 
that 5 — 10% of satellite systems are as flat as the MW’s 
11 classical satellites. When it is required that the velocities 
of at least 8 of the 11 satellites should point within the 
narrow angle claimed by Pawlowski & Kroupa (2013) for the 
MW satellites, this fraction decreases to ~0.1% (Pawlowski 
et al. 2014). In the case of the M31 thin satellite plane, 
Bahl & Baumgardt (2014) found that, while similar spatial 
distributions of satellites are quite common in ACDM, there 
is only a 2% chance that 13 out of the 15 members in the 
plane would share the same sense of rotation. In similar 
studies, Ibata et al. (2014c) and Pawlowski et al. (2014) 
found an even lower occurrence for the M31 plane, with only 
~0.1% of ACDM systems having a similar configuration. 

Extending the above analysis to galaxies outside the LG 
is constrained by observational limitations, but some addi¬ 
tional tests can be performed. Cautun et al. (2015) stud¬ 
ied the flattening of the satellite distribution around iso¬ 
lated central galaxies in the SDSS, as viewed on the plane 
of the sky, finding good agreement between data and cos¬ 
mological simulations. Using a similar approach, Ibata et al. 
(2014b) claimed a higher degree of flattening in their data, 
but their conclusions may be affected by systematics (see 
Cautun et al. 2015, for a discussion of this study). The ex¬ 
pected signature of planar rotation has been investigated 
by considering the velocity correlation of satellite pairs ob¬ 
served on opposite sides of the host galaxy. Initially, using a 
sample of 23 systems, Ibata et al. (2014a) reported a signifi¬ 
cant excess, when compared to ACDM predictions. Cautun 
et al. (2015, see also Phillips et al. 2015) found that this ex¬ 
cess decreases rapidly as the sample size is increased and also 
that the expected mirror image signal is absent for satellite 
pairs on the same side of the host galaxy, suggesting that 
the claimed excess is not robust. 

In this study, we examine planar configurations of satel¬ 
lites identified in MW- and M31-like mock ACDM cata¬ 
logues and compare them to the planar structures observed 
in the LG. Among others, we revisit the claims by Ibata et al. 
(2014c) and Pawlowski et al. (2014) that the two planes of 
satellites found in the LG are inconsistent with AGDM pre¬ 
dictions. Along with similar studies, those works are based 
on two important axioms: that the majority of planar config¬ 
urations are the same, and that the planes detected around 
the MW and M31 are representative of planar structures in 
general. We will examine these two conjectures within the 
context of ACDM predictions. Starting from high resolution 
cosmological simulations, we create mock catalogues that 
are used to identify planar satellite configurations similar to 
those found in the LG. We proceed to study the properties 
of the most prominent planes of satellites and to compare 
those ACDM predictions with the two planes of satellites 
observed in MW and M31. 

This paper is organized as follows. In Sec. 2 we intro¬ 
duce the cosmological simulations as well as the selection 
criteria used to identify MW and M31 analogue systems; in 
Sec. 3 we present an objective method to identify spatially 
and kinematically rare planes; in Sec. 4 we identify planes 
of satellites for PAndAS-like mocks; in Sec. 5 we analyse 
MW-like systems; we conclude with a short discussion and 
summary in Sec. 6. 


2 DATA AND SAMPLE SELECTION 

This work makes use of two high-resolution ACDM DM- 
only cosmological simulations: the Millennium-II (MS-II; 
Boylan-Kolchin et al. 2009) and the Copernicus Complexio 
(COCO; Hellwing et al. 2015) simulations. Instead of us¬ 
ing the original MS-II, which was run assuming a Wilkin¬ 
son Microwave Anisotropy Probe (WMAP)-l cosmogony, 
we employ a modified version of the simulation that has 
been rescaled to the WMAP-7 cosmology (flm = 0.272, 
Ha = 0.728, (Tg = 0.81 and n* = 0.968) using the scaling al¬ 
gorithm of Angulo & White (2010, see also Guo et al. 2013). 
The rescaled MS-II corresponds to a simulation in a periodic 
box of side-length 104.3 /i“^Mpc containing 2160^ DM par¬ 
ticles, each particle having a mass, nip = 8.50 x 10® /i“^M©. 

COCO simulates a smaller, roughly spherical, cosmo¬ 
logical volume, V = 2.25 x 10^ (/i“^Mpc)®, equivalent to a 
cubic box of side-length 28.2 /i“^Mpc, but at much higher 
resolution than MS-II, having 2374® particles each with 
mass, nrip = 1.14 x 10® /i“®M©, and a gravitational soft¬ 
ening length, e = 0.23 /i“^kpc. This volume is embedded 
within a larger periodic box, of side-length 70.4 /i“^Mpc, 
that is resolved at a significantly lower resolution (see Hell¬ 
wing et al. 2015, for more details). COCO uses the same 
WMAP-7 cosmological parameters as the rescaled MS-II. 

To construct mock catalogues, we use the semi-analytic 
galaxy formation model of Guo et al. (2011) that has been 
implemented for the rescaled version of MS-II in Guo et al. 
(2013) and for COCO in Guo et al. (2015). The semi- 
analytic model has been calibrated to reproduce the stellar 
mass, luminosity and autocorrelation functions of low red- 
shift galaxies as inferred from SDSS. The abundance and ra¬ 
dial distribution of satellite galaxies predicted by the model 
are in good agreement with SDSS data (Wang & White 2012; 
Wang et al. 2014), which makes the Guo et al. model a good 
test bed for studying planar structures of satellites. 

Due to the relatively low resolution of MS-II, many of 
the satellite galaxies of interest for this work are found in 
haloes close to the resolution limit of ~10® /i“®M©, which 
raises questions about the accuracy of the properties and 
orbits of these objects, especially after infall into the main 
halo. To test for any potential systematic effects arising from 
the limited resolution of MS-II, we compare with the results 
of COCO. Any such effects are significantly reduced or even 
absent in COCO, since it has 75 times higher mass resolu¬ 
tion and four times better spatial resolution. 

We select counterparts to the two massive members 
of the Local Group by identifying DM haloes with similar 
masses in the range (1 — 3) x 10®®M©, which is consistent 
with the mass of the MW and M31 halo (e.g. Fardal et al. 
2013; PifH et al. 2014; Cautun et al. 2014b,a; Veljanoski et al. 
2014; Gonzalez, Kravtsov & Gnedin 2014, for a compilation 
of other measurements and discussions of systematic effects 
see Courteau et al. 2014; Wang et al. 2015). This results in 
2849 MS-II haloes and in 63 COCO haloes in the required 
mass range. Compared to the previous studies that anal¬ 
ysed the incidence of the M31 plane of satellites (e.g. Bahl 
& Baumgardt 2014; Pawlowski et al. 2014), we adopted a 
broader mass range to account for the large uncertainty in 
the mass measurements and also for possible systematic ef¬ 
fects. We checked that the exact mass range used does not 
affect our final results. 
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Figure 1. Illustration of the procedure to identify planes of satellites. The panels show the same system, with the small open symbols 
corresponding to satellites identified around a central galaxy (large symbol in the centre). We take every subsample of any given number 
of satellites and compute the plane corresponding to that configuration. The panels shows three candidate configurations with Afgat = 7, 
11 and 15 members (red filled symbols) out of a maximum of A^max = 20 satellites. The dashed line corresponds to the best-fitting plane 
and the two dotted lines show the thickness of the plane, rj_ (see Eq. 1). To determine which configuration stands out the most, we 
compute the plane prominence, IP. It specifies that the plane has a probability of one in IP to be a statistical fluctuation, e.g. plane 2 
corresponds to one chance in 660 to be caused by a fluctuation. Out of the three, the rarest or most prominent plane is number 2 since 
it has the largest prominence. 


3 IDENTIFYING PLANAR 

CONFIGURATIONS OF SATELLITES 

In this section we introduce an objective method to iden¬ 
tify, for each halo, the rarest plane of satellites, both spa¬ 
tially and kinematically. The method works by identifying 
the subsample of satellites that is the least likely to be ob¬ 
tained by chance. This is motivated by recent observations 
that have shown that only a subset of the satellite galaxies 
are potentially distributed along a plane. For example, out 
of the 27 M31 satellites in the PAndAS footprint, a signifi¬ 
cant plane is found for 15 of them (IbatalS), while the entire 
population is no more planar than would be expected for an 
isotropic distribution of equal size (Conn et al. 2013). Tully 
et al. (2015) found that the Centaurus A Group shows evi¬ 
dence for two planes of galaxies, which, between them, con¬ 
tain 27 out of the 29 members with known distances. Even 
for the MW, while the 11 classical satellites are found on a 
plane, only 8 have orbital poles in a narrow angle indicating 
a possibly long lived planar structure (Pawlowski & Kroupa 
2013). 

Previous methods of identifying planar configurations 
have been based on first examining the data and only in a 
second step defining an approach for selecting planes, result¬ 
ing in a selection method that is both subjective and a pos¬ 
teriori. For example, the 15 members of the M31 plane have 
been found by noticing that, when increasing the number 
of satellites associated with a plane, the best-fitting plane 
hardly changes and that the thickness of the plane increases 
only slowly. If one considers more than 15 satellites, this then 
leads to a larger change in the thickness of the best-fitting 
plane and in its direction (Ibatal3). Applying such a method 
to a large sample of systems is undesirable, since it implies 
choosing at least two threshold parameters: the maximum 
allowed changes in the thickness of a plane and in its direc¬ 
tion when adding an extra valid plane member. There are no 
a priori well-motivated values for those threshold parame¬ 
ters. In addition, the two thresholds should likely depend on 
the radial distribution of satellites, since more radially con¬ 


centrated distributions will likely have thinner planes. While 
other methods of identifying planar distributions have been 
proposed (e.g. Conn et al. 2013; Gillet et al. 2015), they all 
involve one or more subjective parameters. In contrast, the 
method we present here does not involve any such param¬ 
eters and, in addition, it naturally takes into account the 
radial distribution of satellites in each system. 


3.1 Spatial planes 


We start by introducing a method for identifying the most 
prominent spatial plane. When comparing planes that con¬ 
tain the same number of satellites, the one that stands out 
the most is the thinnest plane. Difficulties arise when we 
have to compare planes with two different numbers of mem¬ 
bers, since it is not trivial to decide which one stands out 
more. In a nutshell, we identify all possible subsamples of 
satellites, out of a total sample of A^ax satellites, and, of 
those, we select the one configuration that is the least likely 
to be caused by statistical fluctuations. This is illustrated in 
Fig. 1. 

In a first step, we identify the satellite subsets that are of 
interest for our study. We do so using the approach detailed 
in Appendix Al. For each such subsample, we find the best¬ 
fitting plane, which is the plane that minimizes the root- 
mean-square distance of the satellites from it. For this, we 
define the plane thickness, r±, as 


ri. 


E£r (npia 


Xi) 




( 1 ) 


where Agat is the number of satellites in the subsample and 
Hpiane deuotes the normal to the plane. With xi we denote 
the position of each satellite in a coordinate system whose 
origin is the central host galaxy. The plane thickness, r±, is 
in fact the dispersion in the distance of the satellites from a 
plane that goes through the central galaxy. The best-fitting 
plane is the one that minimizes r_L. The normal to this plane. 
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Hpiane, IS given by the eigenvector corresponding to the low¬ 
est eigenvalne of the inertia tensor of its members. 

Each resulting plane is characterized in terms of its 
prominence, 7, such that, the larger the prominence, the 
least likely it is that the plane is due to a chance alignment. 
For example, for plane i that has Algat; i members and a 
thickness, r±- i, the spatial prominence is defined as: 

rpplane i _ _ 1 _ 

spatial I iVsat;i) 

where the denominator gives the probability of obtaining 
by chance a configuration of A^sat; i satellites that is thinner 
than rj_; i. This probability is computed using 10® isotropic 
distributions of satellites as outlined in Appendix A3. Since 
the radial distribution of satellites has a strong effect on the 
thickness of the resulting planes, we generate each isotropic 
realization to have the exact same radial distribution as the 
system under study. 

Now it is only natural to characterize the most promi¬ 
nent, or rarest, plane as the one that is the least likely to be 
obtained by chance. Using our notations, this can be written 
formally as: 

m rarest mplane i /o\ 

IPspatial = max , (3) 

all planes i L J 

which says that the rarest spatial plane is the one that has 
the largest prominence. It is important to note that, within 
this approach, every halo contains a rarest plane. Determin¬ 
ing if this rarest plane is statistically significant is a separate 
question that we will address in Sec. 4. 


3.2 Spatial and 2D-kinematic planes 

The observational data for the MW and M31 satellites con¬ 
tains both positions and velocity information for these ob¬ 
jects. It is natural to try to incorporate this additional veloc¬ 
ity information into the detection of planar configurations of 
satellites. The M31 satellites have only radial velocity mea¬ 
surements, so the full 3D velocities are unknown. But since 
the M31 plane of satellites is almost parallel to the line of 
sight, the radial velocities can be used to estimate the sense 
of rotation of each satellite with respect to the best-fitting 
plane. In the following, we describe how to select spatial + 
2D-kinematic planes, which are at the same time spatially 
thin and have a large number of members that share the 
same sense of rotation. 

Before continuing, it is important to discuss some po¬ 
tentially misleading nomenclature used in previous studies. 
Satellites sharing the same sense of rotation have been re¬ 
ferred to as corotating satellites (e.g. Ibatal3, Bahl & Baum- 
gardt 2014). This nomenclature is confusing since corotation 
is normally used to denote a rotation around a common axis. 
Thus, two satellites corotate if their orbital poles are very 
close together. In the absence of 3D velocities, we only know 
that, when projected on the best-fitting plane, 13 out of the 
15 satellites rotate in the same sense, either clockwise or 
counter-clockwise. 

In addition to the steps described in Sec. 3.1, for each 
satellite subset we also determine the number of members 
that share the same sense of rotation relative to the best¬ 
fitting plane. To determine the direction of rotation of each 


member, we take the scalar product between the plane nor¬ 
mal and the orbital momentum of the satellite. A posi¬ 
tive scalar product corresponds to clockwise rotation, while 
a negative one corresponds to counter-clockwise rotation. 
The number of satellites sharing the same sense of rotation, 
Ns.s.i., is the maximum of the number of objects rotating 
clockwise and those rotating counter-clockwise. Following 
this step, we assign to each plane a 2D-kinematic promi¬ 
nence, CP 2 D-kin, which is defined as: 


rpplane i 

‘^2D-kin 


1 

PO Na.s.v.; i I A^sat; i) 


( 4 ) 


where the denominator gives the probability of obtaining 
by chance a configuration of Wat; i satellites in which at 
least Ns.s.i.-, i of them share the same sense of rotation. The 
procedure for computing this probability is detailed in Ap¬ 
pendix A3. 

We define the rarest spatial + 2D-kinematic planes as 
the one whose spatial and 2D kinematical distribution are 
the least consistent with a statistical fluctuation. Thus, 


rp rarest 

spatial + 2D-kin 


max 

all planes i 


rpplane i rpplane i 
'^spatial ^ *^20 —kin ’ 


( 5 ) 


that is the plane that maximizes the product of the spatial 
and the 2D-kinematic prominences. 


3.3 Spatial and 3D-kinematic planes 


In the case of the MW, the 3D velocities of the 11 classical 
satellites are known. This suggests that for the MW sys¬ 
tem one can identify planes that are both spatially thin and 
show a large degree of coherent 3D kinematics. For this, one 
needs to construct a cost function that rewards systems in 
which most satellites have orbital poles close together and 
penalizes systems in which the orbital poles are isotropically 
distributed. For example, to study long lived planar config¬ 
urations, the cost function would preferentially reward sys¬ 
tems in which the orbital momentum of its members is close 
to parallel or anti-parallel to the normal to the best-fitting 
plane. For this work, we employ a variant of the cost function 
suggested by Pawlowski & Kroupa (2013), since this one has 
been used in other studies that claim a tension between the 
MW satellite plane and ACDM predictions (e.g. Pawlowski 
et al. 2014). That function has been proposed after examin¬ 
ing the orbital data of the MW satellites, and as such it is an 
a posteriori definition of the rotation characteristics of the 
Galactic satellite distribution. It may be that other satellite 
planes in the Universe have different orbital structures, in 
which case that cost function may not be optimally suited 
for characterizing their kinematical structure. 

To compute the amount of kinematical information, we 
proceed as follows. For each of the satellite subsets used in 
Sec. 3.1, we determine the dispersion in the orbital poles, 
i.e. directions of the orbital momenta, of its members as 


Astd — 


arccos^ (Uorbit • nprbit; i 

Wat 


( 6 ) 


where Uorbit; i denotes the orbital momentum direction of 
each member of the plane and Horbit the mean orbital pole 
of all the Wat members found in the plane. Compared to our 
approach, Pawlowski et al. (2014) applied Eq. (6) to only 8 
out of the 11 satellites found in the MW satellite plane. This 
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Table 1. The characteristics of the rarest planar configuration of satellites around the M31 and the MW. 


ID 

Plane type 

-Vsat 

Ns.s.t. 

r± 

(kpc) 

Til 

(kpc) 

Astd 

(degrees) 

y rarest 

fraction of systems with 
more prominent planes (%) 









ACDM 

isotropic 

M31 plane of satellites 

(1) 

spatial 

14 

11“ 

10 

220 ^ 2 ! 

- 

X 10^ 

12l^ 

1 2"^^'^ 

-*^•^-0.6 

(2) 

spatial -h 2D-kinematic 

15 

12“ 

12 5+0 '^ 

214+^^ 

^-*^^-23 

- 

3.4^^;® X lO"* 


0 24+0-^® 
'^•'^^-0.14 

MW plane of satellites 

(3) 

spatial 

11 

8'’ 

20 7+*^-^ 

• _o.6 

129^4 

63l® 

X 10^ 

I 2 I) 

2 

^•^-0.4 

(4) 

spatial -h 2D-kinematic 

11 

8 ’’ 

• _o.6 

129^4 

63l® 


IOI 5 

2.oli;f 

(5) 

spatial -h 3D-kinematic 

11 

8 ^ 

1 _QQ 

129t^ 

63l| 

8.311^0 X 10® 

S-Olay 

^•^®-0.39 


The M31 plane was found using 27 satellites within the PAndAS footprint. The MW plane was found using the brightest 11 classical 
satellites and a Galactic obscuration zone of 33% of the sky. The table columns give: a plane ID for easy reference, plane selection 
method, the number of satellites in the plane (A^sat) and how many of those share the same sense of rotation with respect to the 
best-fitting plane(As.s.r.)> the plane thickness (rj^) and radial extent (t||), the orbital pole dispersion for the plane members (Agtd)? 
the prominence of the plane (y rarestfraction of systems with more prominent planes, i.e. higher Trarest^ £qj. ^.nd isotropic 

satellite distributions. 

“ The sense of rotation of And XXVI is highly uncertain since it depends on the object’s radial distance from the MW. And XXVI 
is more likely to counter-rotate (55% probability) than to rotate in the same sense as the majority of the satellites in the plane (45% 
probability). This is in contradiction with the Ibatal3 results, who claimed that And XXVI is one of the 13 plane members that share 
the same sense of rotation. We suspect that Ibatal3 calculated the sense of rotation using the distance corresponding to the peak of the 
radial distance PDF, which indeed would result in the claimed result. The radial distance PDF of And XXVI is highly asymmetrical, so 
the position of the PDF peak does not necessarily characterize the most likely outcome. 

^ Due to large proper motion errors, there are uncertainties in the sense of rotation of: Sextans, Carina, Leo I and Leo 11. This mean 
that the number of satellites sharing the same sense of rotation is 6, 7, 8 or 9 with a probability of 1, 16, 50 and 33%, respectively. 


was motivated by the observation that only 8 out of the 11 
members show close orbital poles^. In principle, we could 
follow a similar approach and take a subset of the plane 
members that shows the most concentrated orbital poles. 
This would amount to taking a subset of a subset, since our 
planes are already subsets of satellites from a larger sample. 
We prefer not to do so since it would add an additional layer 
of complexity to this method and also a significantly higher 
computational cost. 

After applying Eq. (6) to each plane, we define the 3D- 
kinematic prominence, IPsD-kin, of plane i as 


rpplane i _ _1_ 

3D — kin ^ \ iTVr V/ 

P (^ Astd; i I Asat; i j 

where Astd; i denotes the orbital dispersion of the plane. The 
denominator gives the probability of obtaining by chance a 
confignration of Ns.ax- i that have an orbital pole dispersion 
less than Aatd; i- The procedure for computing this proba¬ 
bility is given in Appendix A3. 

Now we can define the prominence of the rarest spatial 
+ 3D-kinematic plane as 


rn rarest 

•^spatial -|- 3D-kin 


max 

all planes 




rpplane i rpplane 

spatial 3D —kin 


in ’ 


( 8 ) 


which is the plane of satellites whose spatial and 3D kine- 


^ Incorporating such a posteriori considerations incurs the danger 
of designing tests that are specifically matched to one particular 
system and that may not be characteristic of the larger popula¬ 
tion. 


matic distribution is the least likely to be a statistical fluc¬ 
tuation. 


4 M31-LIKE PLANES OF SATELLITES 

Here we investigate the characteristics of the rarest planes of 
satellite galaxies as found in mock PAndAS-like catalogues. 
The goal is to obtain a better understanding of the M31 
plane of satellites and to compare it to the ACDM predic¬ 
tions. 

To create PAndAS-like mocks, we use the host halo sam¬ 
ple described in Sec. 2. For each of those hosts, we start 
by finding all the satellites with stellar masses larger than 
2.8 X lO'^M© (as proposed by Bahl & Baumgardt 2014) that 
are within a radial distance of up to 500 kpc. To reproduce 
the observational geometry, we place the observer at a dis¬ 
tance of 780 kpc from the centre of the host halo, which is 
equivalent to the MW-M31 distance (Conn et al. 2012). For 
each satellite identified earlier, we compute its sky coordi¬ 
nates, as seen by the observer. Out of all the satellites within 
the PAndAS mask that are also more than 2.5° from the host 
(Ibatal3), we keep only the 27 objects that have the largest 
stellar masses. If there are fewer than 27 satellites within the 
required geometry, we discard that host halo. For each host, 
we place the line-of-sight along three mutually perpendicu¬ 
lar directions consisting of the simulation’s x, y and z axes. 
Due to the highly asymmetrical PAndAS volume, this will 
result in somewhat different satellite distributions, hence in¬ 
creasing the overall statistics. After applying this procedure. 
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we end up with 7350 mock satellite systems in MS-II and 
180 in COCO. 


4.1 The M31 system 

In a first step, we apply our method to the actual PAndAS 
observations of M31. We do not use the same plane identifi¬ 
cation method as Ibatal3, so it is important to check what it 
is that our approach identifies as the most prominent plane 
of the M31 system. To account for observational errors, we 
generate 1000 Monte Carlo realizations that sample the ra¬ 
dial distance PDFs (Conn et al. 2012, Table 1) and radial 
velocities (Collins et al. 2013, Table 5) of the M31 satellites. 

In the case of the M31 system, the observational data 
allow for the identification of spatial and spatial + 2D- 
kinematic planes. For each Monte Carlo realization of the 
M31 system, we identify the rarest plane. Due to the large 
radial distance errors, the rarest plane can vary between re¬ 
alizations. For example, the rarest spatial plane contains 14 
members in 72% of the cases and 13 members in 22% of 
realizations. In the remaining 6% of realizations it contains 
even fewer satellites. For simplicity, we take the rarest plane 
as the one that is identihed as such in the largest number 
of realizations. The rarest planar configurations of the M31 
system and its characteristics are shown in Table 1. 

We find that the rarest spatial plane consists of 14 satel¬ 
lites that are the same as the 15 members of the plane re¬ 
ported by Ibatal3, except And III. This is in agreement with 
the results of Ibatal3, who point out that choosing 13 or 14 
satellites results in a higher spatial significance, i.e. lower 
probability of being obtained from an isotropic distribution, 
than for the full sample of 15. The spatial + 2D-kinematic 
plane found by our approach is the same as the one reported 
by IbatalS, even though our plane identification method is 
different. IbatalS reported the significance of the M31 plane 
as compared to an isotropic distribution, so it is possible that 
they inadvertently choose the parameters of their method 
(see Sec. 3) such that it maximizes the plane significance. 
If that was the case, then both plane finding methods are 
basically the same. 


4.2 The rarest M31-like planes 

To better understand the M31 plane of satellites, we start 
by assessing the chance of obtaining more prominent planar 
configurations within ACDM. This is shown in Fig. 2, where 
we plot the cumulative distribution function (CDF) of the 
prominence, IPa’^atfaii of fh® rarest spatial planes. There is 
a very good match between the MS-II and COCO haloes, 
which suggests that satellite planes found in MS-II are not 
signihcantly affected by the limited resolution of the sim¬ 
ulation. We find that most haloes have prominent planar 
configurations, for example 37, 12 and 4% of haloes have 
planes with ^ 10^ i 10^ s-od 10^ respectively. 

The prominence is not simply the inverse of the proba¬ 
bility for isotropic satellite distributions for reasons that will 
be discussed in Sec. 4.3. As such, the figure also shows the 
rarest planes found in an isotropic distribution of satellites. 
To obtain these, for each of the ACDM haloes, we gener¬ 
ate an isotropic distribution by choosing random polar and 
azimuthal angles, while keeping the radial position of each 


significance [o] 



Figure 2. The CDF of the prominence, of the rarest 

spatial plane of satellites for mock PAndAS observations. The 
solid line gives the MS-II results, while the dashed line shows re¬ 
sults for COCO, which has much higher resolution. The dotted 
curve gives the expectation for isotropic satellite distributions. 
The vertical dashed line and shaded region show the prominence 
and the Icr error for M31’s spatial plane of satellites. We find that 
12tf% of ACDM haloes have a more prominent planar configu¬ 
ration than M31. The top x-axis shows the detection significance 
of each plane computed using the isotropic CDF (dotted curve) 
and accounts for the look-elsewhere effect. 


significance [o] 



Figure 3. As Fig. 2, but for the prominence, T, 
of the rarest spatial + 2D-kinematic plane of satellites. This case 
corresponds to the M31 plane identified by Ibatal3. For this test, 
8.8^^ g% of ACDM hosts have a more prominent plane than M31. 

satellite fixed. We then apply the same plane identification 
procedure to each isotropic satellite distribution. Unsurpris¬ 
ingly, we find a clear difference between the isotropic and 
ACDM results, with the isotropic CDF shifted towards the 
left of the ACDM CDF. This suggest that, compared to a 
uniform distribution, there is more structure in the distribu¬ 
tion of ACDM satellites, in agreement with previous studies 
(e.g. Libeskind et al. 2005; Wang, Frenk & Cooper 2013; 
Pawlowski et al. 2014). 

The corresponding M31 plane, entry (1) from Table 1, 
is shown as the dashed vertical line. It has a prominence, 
^spatfai = l-OlJ's ^ > which means that for an isotropic 

distribution there is only a 1 in 1000 chance of obtaining 
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a thinner plane with 14 members. This result was found by 
computing the prominence of the rarest plane for each Monte 
Carlo realization of the M31 system. Following this, we quote 
the median value and the Icr interval, corresponding to the 
16th to 84th percentiles. We find that the spatial plane of 
the M31 system is consistent with ACDM expectations, since 
there is a chance of finding an even more prominent 

plane in ACDM. In fact, the M31 spatial plane is consistent, 
at 2.5 (t, even with an isotropic distribution, since there are 
1.2i'iJ'g% more extreme systems in this case (these results 
are summarized in Table 1). 

Fig. 3 shows the prominence of the rarest spatial + 2D- 
kinematic planes. The general conclusions are the same as 
for the previous figure: we find a very good match between 
the MS-II and the COCO results; and ACDM satellite dis¬ 
tributions have more prominent planes than isotropic distri¬ 
butions. The bumpy aspect of the CDF curves for the MS-II 
and the isotropic case reflects the discrete nature of the 2D- 
kinematical test, since the number of plane members sharing 
the same sense of rotation always takes integer values. For 
this test, the corresponding M31 plane, entry (2) from Ta¬ 
ble 1, which is the one identified by Ibatal3, is characterized 
by a prominence, T™riai + 2 D-kin = S.dlJ J x 10'*. There are 
8 . 8 ^ 13 % ACDM haloes with more prominent planes sug¬ 
gesting that the M31 plane of satellites is in agreement with 
ACDM predictions. 

4.3 The detection significance of a plane 

We now discuss the detection significance of a plane of satel¬ 
lites, i.e the probability that a plane is due to a statistical 
fluctuation. For this, we need to take into account the “look- 
elsewhere” effect. This is a phenomenon in statistics where 
an apparently statistically significant observation may have 
actually arisen by chance because of the large size of the 
parameter space to be searched. It represents an important 
effect for cases where one does not have an a priori model 
or prediction to where the signal should appear, and, hence, 
when one needs to search for a signal in a large range. In 
such cases, the significance calculation must take into ac¬ 
count that a high statistical fluctuation anywhere in that 
range could also be considered as a signal (e.g. see Gross & 
Vitells 2010, for a more rigorous discussion). The effect is 
particularly relevant in particle physics, and, in general, in 
any field in which one searches for uncommon events. 

The look-elsewhere effect is important since we do not 
know a priori what is the number of satellites that we ex¬ 
pect to hnd in a plane. In the case of PAndAS-like obser¬ 
vations, there are 27 satellites in each system, which means 
that the most prominent plane can have anywhere between 
3 to 27 members. To estimate the significance of a plane, we 
need to compute the probability that a statistical fluctuation 
generating such a prominent plane appeared for any com¬ 
bination of 3,4,..., 27 satellites — this is the marginalized 
probability. This is different from the conditional probabil¬ 
ity that has the number of satellites chosen as a prior and 
whose inverse gives the prominence of a plane (e.g. Eq. 2). 
The marginalized probability is the isotropic CDF shown in 
Figs. 2 and 3. In other words, the probability that a plane is 
a statistical fluctuation is given by the fraction of isotropic 
realizations that have a more prominent plane. Using this, 
we compute the significance of each plane, which is shown as 



Figure 4. The trial factor for the look-elsewhere effect as 
a function of the prominence of the rarest plane. The ratio, 
trial factor/?gives the marginalized probability of obtain¬ 
ing by chance (i.e. in an isotropic distribution) a plane with promi¬ 
nence, yra.rest regions of interest are roughly ?/patfai ^ 

and if/patfai + 2D-kin ^ Corresponding to the intervals in 

Figs. 2 and 3 where the CDF is different from unity. 

the top x-axis of Figs. 2 and 3, with the tick marks spaced 
at 0.2 (t intervals. The significance is expressed in multiples 
of the standard deviation, cr, of a normal distribution. Note 
that the one-to-one map between the prominence (bottom 
x-axis) and the significance level (top x-axis) differs between 
the two figures. 

Some previous studies have incorrectly referred to the 
conditional probability as the detection significance of planes 
of satellites (e.g. Kroupa, Theis & Body 2005, Ibatal3). To 
better emphasize the difference between the two, we intro¬ 
duce the concept of trial factor (e.g. see Gross & Vitells 
2010). This is the ratio between the marginalized and the 
conditional probability to obtain a statistical fluctuation 
with prominence, T. The former corresponds to the actual 
detection significance and is given by the GDF of an isotropic 
distribution. The latter corresponds to since this is the 
definition of a plane’s prominence (e.g. Eq. 2). Thus, 

trial factor - y rarest^ 

Liitxi IdLLUl — -T- 

rarest^ 

= GDEi.otropic(^ . (9) 

For example, the M31 spatial -h 2D-kinematic plane has a 
trial factor of 115 (see entry (2) in Table 1 for numerical 
values). Thus, the chance of it being a statistical fluctuation 
is 115 times higher than naively expected if one considers 
only random planes with 15 members. Inevitably, this means 
that Ibatal3 have overestimated the detection significance 
of the M31 plane by more than two orders of magnitude. 
For an isotropic distribution there is a 0.34% probability 
of obtaining a more prominent plane, and, hence, the M31 
plane corresponds to a 2.9 (j detection. 

In Fig. 4 we show the trial factors for spatial and spatial 
-h 2D-kinematic planes. The regions of interest are roughly 
3’a''patTai > 10 and TapItrai + 2D-kin ^ 100 Corresponding to the 
intervals in Figs. 2 and 3 where the isotropic GDF is differ¬ 
ent from unity. In those intervals, the trial factors increase 
only slowly with the plane prominence, so, to a first approx- 
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imation, the two plane types have a trial factor of ~10 and 
~100, respectively. The spatial + 2D-kinematic planes have 
a higher trial factor due to the larger range used to search 
for the most prominent plane, since, on top of the spatial 
distribution, also the 2D kinematics are considered. 


4.4 The characteristics of rare planes 

In Fig. 5 we show the characteristics of the rarest spatial + 
2D-kinematic planes that are at least as prominent as the 
M31 plane. We describe the planes in terms of the number 
of members sharing the same sense of rotation, the plane 
thickness and the radial extent of the plane, ry. This latter 
property characterizes the dispersion of the satellites within 
the plane and it is calculated as the mean sum of the squares 
of the distance projected on to the best-fitting plane. We 
choose these plane characteristics to be consistent with pre¬ 
vious studies that investigated the incidence of the MW and 
M31 plane of satellites in term of these properties (e.g. Bahl 
& Baumgardt 2014; Ibata et al. 2014c). 

Fig. 5 shows that there is considerable variation among 
the properties of the most prominent planes, suggesting that 
each plane is different. For example, the number of members 
sharing the same sense of rotation, Ns.s.r., can take values be¬ 
tween Wat/2 to Wat. The top panel of the figure shows that 
the planes have W.s.r. values spanning the full allowed range, 
although there is a higher preference for W.s.r. — Wat, since 
that will result in a higher prominence. The behaviour in the 
middle panel is governed by two requirements. Firstly, to be 
prominent, planes with a small number of members need to 
be very thin since such structures cannot have a high 2D- 
kinematic prominence, which explains the distribution seen 
in the left-half of the panel. Secondly, once the number of 
members is high enough, ~16 in this case, the 2D-kinematic 
prominence can be by itself very large, so that such planes 
do not necessarily need to be very thin. This explains the 
large scatter in the r± values seen in the right-half of the 
middle panel. And lastly, the prominence of a plane does not 
depend on ry, which explains the large scatter in ry values 
seen in the bottom panel of the figure. 

Fig. 5 also shows the properties of the M31 plane of 
satellites whose position is marked with a large triangle. 
The M31 plane is within the scatter expected for ACDM 
planes, although it does stand out as having an unusually 
large radial extent, ry. 


4.5 The incidence of rare planes 

For each halo, we study the incidence of the rarest spatial + 
2D-kinematic plane among the distribution of satellites of 
all other ACDM haloes. This is motivated by the studies of 
Ibata et al. (2014c) and Pawlowski et al. (2014) that inter¬ 
preted the low incidence of the M31 plane as evidence for 
an inconsistency between observed planes of satellites and 
ACDM predictions. 

We define the incidence or frequency of a plane using the 
approach of Ibata et al. (2014c). Each plane of satellites is 
characterized by; its number of members, Wat; how many of 
them share the same sense of rotation, W.s.r.; the thickness, 
r^; and radial extent, ry, of the plane. Then, the frequency 
or incidence, /ACDM(Wat, ^ry, ^W.s.r.), of this plane 
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Figure 5. The characteristics of the rarest spatial -h 2D- 
kinematic planes of satellites that are at least as prominent as 
the M31 plane, i.e. 2 D-kin ^ 1*^ ^ 1*^^‘ There are 852 

such systems in MS-II. The grid cells are coloured according to 
the number of systems with those properties. The three panels 
show the number of satellites sharing the same sense of rotation, 
Ns.s.t. (top)? the thickness, rj_ (centre), and the radial extent, 
ry (bottom), of the best-fitting plane, as a function of the num¬ 
ber of satellites in the plane, A^sat- The large triangle shows the 
properties of the M31 plane. 



is given by the fraction of ACDM systems that have a similar 
plane. To describe the procedure, we exemplify it for the case 
of two systems A and B. We are interested in the frequency 
of the rarest plane of satellites of halo A and we wish to find 
out if halo B has a similar plane. We take all possible satellite 
configurations of system B that have members, which 
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Table 2. The incidence of the M31 and the MW plane of satellite galaxies. 


ID 

Host 

/ACDM(.^^sat, ^rj^, ^As.a.r.) 

/ACDM(.^aat, , sgAgtd) 

fraction of ACDM systems 
with lower frequencies {%) 

(2) 

M31 

e.sta.r ^ 10-4 

- 


(4) 

MW 


- 


(5) 

MW 

- 

l.sli;® X 10-3 

to b- 

+ 1 


The table columns give: the plane ID from Table 1, the central galaxy (M31 or MW), the incidence, /ACDM(^sat, JSr||, ^A^a.s.r.) 
and /ACDM(^aat, ^r||, ^Astd), of similar planes among the MS-II haloes, and the fraction of ACDM haloes that have planes of 
satellites with even lower frequencies. 



Figure 6. The CDF of the incidence, /ACDM(fVsat, ^rx, ^r||, 
^fVs.s.r.), of similar planes to the rarest one in ACDM. For each 
ACDM halo, we take the most prominent plane and find its fre¬ 
quency among all other ACDM haloes (see the text for more 
details). Most planes have a low incidence (half of systems have a 
frequency of 0.02 or lower) indicating that each ACDM halo has 
a different planar configuration. The vertical dashed line and the 
grey area show the incidence and the la error for the M31 plane. 
We find that of ACDM systems have a lower frequency 

than the M31 plane. Thus, the low incidence of the M31 plane is 
not in tension with ACDM, as claimed by Pawlowski et al. (2014), 
but instead is consistent with ACDM expectations. 


we find using the procedure described in Appendix Al. If any 
of those configurations is similar to plane of A, i.e. r5 ^ rf 
and r® rjj^ and ^ ^ , then halo B has at least 

one planar configuration similar to that of system A. We 
compute the frequency for each MS-II and COCO halo, 
by taking the characteristics of the rarest plane identified 
around each halo. Each such rarest plane is compared to the 
satellites distribution of all the other PAndAS-like mocks. 
This equates to comparing the rarest plane found in one halo 
with all possible planes, not only the rarest ones, around all 
other systems. 

The frequency of the rarest spatial + 2D-kinematic 
planes of satellites is shown in Fig. 6, where we plot the CDF 
of the frequency for MS-II and COCO planes. We find that 
over half of ACDM haloes have a plane with a frequency 
of 0.02 or lower, and one tenth have a frequency as low as 
0.001. It is important to note that the low frequency is not 
a consequence of our plane identification method. For each 
halo, our method selects the plane that is the least likely to 
be a statistical fluctuations, which is fully independent of 



Figure 7. The correlation between a plane’s promi¬ 
nence, 3’spatial-|-2D-kin’ i*® incidence, /ACDM(A^sat, ^rx, ^r-ji, 
ISA^s.s.r.), among other ACDM haloes. The small dots correspond 
to the rarest plane found for each MS-II halo. The solid line shows 
the median trend. The large symbol with error bars shows the po¬ 
sition of the M31 plane. The small down-pointing arrows found 
at the bottom of the graph show upper limits corresponding to 
planes that do not have another counterpart in MS-II. 

the planes found around other ACDM haloes. Thus, there is 
a large diversity of planes of satellites. In other words, if we 
find a planar configuration around one system, it does not 
tell us anything about the properties of planes around other 
haloes. 

We performed the same calculation for the M31 system. 
For each Monte Carlo realization, we compute the incidence 
of the rarest spatial + 2D-kinematic plane of that realization 
by comparing with the PAndAS-like mocks. We found that 
the M31 plane has a frequency of 10“^ (Icr confi¬ 

dence interval; see Table 2), which is in good agreement with 
the results of previous studies (Ibata et al. 2014c; Pawlowski 
et al. 2014). This low incidence of the M31 plane has been 
claimed by Ibata et al. (2014c) and Pawlowski et al. (2014) to 
be a source of discrepancy with ACDM. From Fig. 6, which 
shows the M31 frequency as a vertical dashed line, we find 
that, within ACDM, 5.llQ g% of systems have planes with 
even lower frequencies. Thus, the low incidence of the M31 
plane does not pose a challenge to the current paradigm, in 
fact, it is consistent with ACDM predictions. 

Fig. 7 investigates which planes are the ones with the 
lowest incidence. For this, we plot the incidence of each 
plane as a function of its prominence, and find a strong 
anti-correlation between the two, albeit with a large scatter. 
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Table 3. The positions and velocities of the 11 classical MW satellites with respect to the Galactic Center. 


Name 

a/(kpc) 

y(kpc) 

.^(kpc) 

14 (km/s) 

Uj;(km/s) 

14 (km/s) 

Sagittarius 

17.1 ± 1.9 

2.5 ±0.2 

-6.5 ±0.5 

234 ± 7 

19 ±21 

224 ±21 

LMC 

-0.5 ± 0.3 

-41.8 ± 1.6 

-27.5 ± 1.1 

-42 ± 12 

-226 ± 13 

234 ± 16 

SMC 

16.5 ± 1.6 

-38.5 ± 2.4 

-44.7 ± 2.8 

2 ± 18 

-161 ± 26 

149 ± 21 

Draco 

-4.4 ±0.3 

62.3 ±4.9 

43.2 ±3.4 

74± 24 

43± 14 

-210 ± 19 

Ursa Minor 

-22.2 ±0.6 

52.0 ±2.1 

53.5 ±2.1 

7± 28 

89 ±20 

-186 ± 20 

Sculptor 

-5.2 ± 0.3 

-9.7 ± 0.7 

-85.3 ±5.9 

-33 ± 44 

188 ± 45 

-99 ±6 

Sextans 

-36.6 ± 1.3 

-56.8 ± 2.6 

57.8 ±2.7 

-168 ± 160 

114 ± 133 

117±127 

Carina 

-25.0 ± 1.0 

-95.8 ± 5.5 

-39.7 ± 2.3 

-74 ±44 

8 ± 19 

40± 41 

Fornax 

-41.1 ± 2.7 

-50.8 ±4.1 

-134 ± 11 

-38 ± 27 

-156 ±42 

113 ± 18 

Leo II 

-77.3 ±4.1 

-58.3 ±3.5 

214 ± 13 

102 ± 127 

237±156 

117 ± 50 

Leo I 

-124 ± 7 

-119±7 

192 ± 12 

-167± 31 

-35 ± 33 

96± 24 


The x-axis points from the Sun towards the Galactic Centre, the y-axis points in the direction of Galactic rotation at Sun’s position and 
the z-axis points towards the North Galactic Pole. Since this is a rotated coordinate system, the uncertainties are correlated and are very 
anisotropic in the plane of the sky. 


The vertical concentrations of points correspond to planes 
that have a very high 2D-kinematic prominence but only 
a very low spatial one, with the discrete nature of the 2D- 
kinematic test leading to many planes having very similar 
+ 2 D-kin values. So, On average, the more prominent 
a plane is, the lower is its incidence among ACDM systems. 
This explains why the M31 plane, which has a high promi¬ 
nence, also has a low incidence. The M31 plane, shown in 
Fig. 7 as a square symbol with error bars, is consistent with 
the object-to-object scatter in the prominence-incidence re¬ 
lation. Compared to the median trend, the M31 plane has 
an ~5 times lower incidence for its prominence, which could 
be due to its unusually large radial extent (see bottom panel 
of Fig. 5). 


5 MW-LIKE PLANES OE SATELLITES 

In this section we investigate the MW system of satellites 
and how it compares with other ACDM planar configura¬ 
tions. Compared to the M31 analysis, there are three main 
differences: we use systems with 11 instead of 27 satellites, 
the survey geometry is different and, most importantly, we 
can perform additional tests since we have proper motion 
data for the MW satellites. In this analysis, we consider 
only the brightest 11 classical Galactic satellites since only 
these objects have measured proper motions. A twelfth satel¬ 
lite, Canes Venitici, has a similar absolute magnitude as the 
faintest of the classical satellites (McConnachie 2012), but 
it does not have a measured proper motion, so we do not 
include it in this study. 

To create MW-like mocks, we consider the 11 satellites 
with the largest stellar masses that are within a distance 
of 260 kpc from the central galaxy (corresponding to Leo 
I which is the furthest at 254 kpc) and that, at the same 
time, are outside an obscuration region consisting of 33% of 
the sky. The obscuration region accounts for the Galactic 
zone of avoidance, where, due to large extinction and confu¬ 
sion by foreground stars, it is possible to have yet undiscov¬ 
ered bright satellites. For this, we use the estimate of Will- 
man et al. (2004) according to which the census of Galac¬ 
tic dwarf galaxies may be 33% incomplete. Yniguez et al. 
(2014) estimates an even higher incompleteness, with most 
of those undetected systems further than 100 kpc from the 


Galactic Gentre. A satellite is lying inside the obscuration 
region, and hence undetected, if its latitude, 6, is in the 
range —Scrit ^ ^ ^ ^'crit, with 0crit = 19.5° (due to an error, 
Wang, Frenk & Gooper (2013) misquoted the critical angle 
as having a value of 9.5°, which would correspond to an ob¬ 
scuration region of 17%). To generate mock observations, we 
take three viewing angles for each simulated halo such that 
the mock north Galactic pole corresponds to the x, y and 
z axes of the simulation. This procedure yields 8547 mock 
satellite systems for MS-II and 189 for COCO. 


5.1 The MW system 

In a first step, we analyse the planar configurations found 
around the MW. For each of the 11 Galactic satellites, we 
take their radial distance and velocity, as well as the er¬ 
rors associated with these quantities, from McGonnachie 
(2012). To compute the mean proper motion, we follow 
the approach of Pawlowski & Kroupa (2013) and weigh 
the different measurements according to their errors. We 
obtain the same mean proper motions as they do, ex¬ 
cept for Draco, where the latest measurement, (/ia,/^^) = 
(0.177 ± 0.063, -0.221 ± 0.063) mas yr"^ (Pryor, Piatek & 
Olszewski 2015), is signihcantly different from the value 
given by Pawlowski & Kroupa (2013). Using this updated 
measurement results in a weighted mean value ~ 

(0.187 ± 0.063, —0.201 ± 0.063) mas yr“^. We also included 
an additional proper motion measurement for Sagittarius of 
(Aa,/i«) = (—2.95 ± 0.21,—1.19 ± 0.16) masyr“^ (Massari 
et al. 2013) that resulted in a mean value of (jlaj'Ps) — 
(—2.711 ± 0.066,—1.043 ± 0.065) masyr“^, nearly the same 
as the mean value used by Pawlowski & Kroupa (2013). 

We transform the satellite positions and velocities to a 
Gartesian coordinate system with the origin at the Galac¬ 
tic Gentre. The x-axis points from the Sun towards the 
Galactic Centre, the y-axis points in the direction of Galac¬ 
tic rotation at the Sun’s position and the z-axis points to¬ 
wards the North Galactic Pole. For this transformation we 
adopt: the distance of the Sun from the Galactic Gentre 
d© = 8.29 ± 0.16 kpc, the circular velocity at the Sun’s 
position, Wire = 239 ± 5 km/s (McMillan 2011), and the 
Sun’s motion with respect to the Local Standard of Rest, 
([/, V, W) = (11.1±0.8,12.2±0.5, 7.3±0.4) km/s (Schonrich, 
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Binney & Dehnen 2010). To account for observational er¬ 
rors, we generate 1000 Monte Carlo realizations of the MW 
system of satellites. We sample the satellite positions and 
proper motions from Gaussian distributions centred on the 
most likely values of each quantity and with dispersion equal 
to the uncertainties. Similarly, we also account for the errors 
in d©, Kirc and the Local Standard of Rest. Following this, 
we transform from heliocentric coordinates to Galactic ones, 
with the result used as input for our plane detection method. 
We summarize in Table 3 the mean positions, velocities and 
Icr errors associated with each Galactic satellite. 

As for the M31 case, for each Monte Carlo realization 
of the MW system we compute the rarest spatial, spatial + 
2D-kinematic and spatial + 3D-kinematic planes. The re¬ 
sults are summarized in Table 1. Independently of the plane 
type, we find that the rarest plane is the one that contains 
all the 11 Galactic satellites. This is in agreement with pre¬ 
vious studies that found that all the classical satellites are 
members of the MW satellite plane (e.g. Kroupa, Theis & 
Body 2005). 

5.2 The rarest MW-like planes 

To put the MW plane of satellites into context, we proceed 
by identifying the rarest planes around mock MW-like ob¬ 
servations. The outcome is shown in Figs. 8-10 that give the 
CDF of the prominence of the rarest spatial, spatial + 2D- 
kinematic and spatial + 3D-kinematic planes, respectively. 
The main conclusions are the same as for the PAndAS mocks 
analysed in Sec. 4.2. 

The look-elsewhere effect again plays an important role 
for the MW-like mocks, with spatial planes having a trial 
factor ~8, while spatial + 2D-kinematic and spatial + 3D- 
kinematic planes have a trial factor ~30. The trial factors are 
roughly constant in the region where the isotropic CDF is 
below unity, reminiscent of the results in Fig. 4. For brevity, 
we do not show these results. The trial factors are lower in 
the case of the MW mocks than in the case of M31, reflecting 
the narrower range used to search for planar conhgurations 
in the MW (11 satellites compared to 27 for the M31). 

Figs. 8-10 also indicate the prominence of the MW plane 
of satellites as a dotted vertical line. The MW plane stands 
out the most in terms of its spatial + 3D-kinematic promi¬ 
nence since in this case it corresponds to a 2.8 ct statisti¬ 
cal fluctuation. Not accounting for the look-elsewhere effect, 
would lead one to estimate the MW plane as a 3.8(J detec¬ 
tion. While the MW spatial + 3D-kinematic plane is quite 
conspicuous, it is consistent with ACDM since 5.0^2']-% of 
galactic-mass haloes have even more prominent planes. 

As we emphasized in Sec. 3.3, one needs to be careful 
when interpreting the results of the 3D-kinematic analysis 
since this test has been designed a posteriori. In fact, one 
could easily come up with other 3D kinematic tests that 
are physically better motivated, as we discussed in Sec. 3.3. 
Given that observationally we have only one satellite sys¬ 
tem with 3D velocities, it is impossible at present to assess 
if the 3D-kinematic test that we have applied is generic and 
thus appropriate to the whole population of satellite sys¬ 
tems or if it is matched to the particular details of the MW 
satellites. If the latter is true, then the fraction of ACDM 
haloes with more prominent planes does not convey any 
physically meaningful information. Thus, there is currently 


significance [o] 



Figure 8. The CDF of the prominence, of the rarest 

spatial plane of satellites for mock MW observations. The vertical 
dashed line and shaded region show the prominence and Icr error 
for the MW plane of satellites, with (12 di 1)% of ACDM haloes 
having a more prominent plane. 


significance [o] 



Figure 9. As Fig. 8, but for the prominence, ^spatfal + 2 D-kin’ 
the rarest spatial -h 2D-kinematic plane of satellites. In this case, 
(10 di 5)% of ACDM hosts have a more prominent plane than the 
MW. 



1 . rnrarest 

plane prominence 

Figure 10. As Fig. 8, but for the prominence, -i- 3 D-kin’ 

the rarest spatial + 3D-kinematic plane of satellites. For this test, 
5.0^2 y9o of ACDM haloes have a more prominent plane than the 
MW.' 
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Figure 11. The characteristics of the rarest planes of satellites 
that are at least as prominent as the MW plane of satellites, i.e. 
^spatfal + 3 D-kin ^ ^ There are 778 such systems. The 

grid cells are coloured according to the number of systems with 
those properties. The three panels show the orbital pole disper¬ 
sion, Astd (top )5 tho thickness, r± (centre), and the radial extent 
of the plane, ry (bottom), as a function of the number of satellites 
in the plane, A^sat- The large triangle shows the corresponding 
characteristics of the MW plane of satellites. 


not enough data to decide if one should be concerned that 
only ~5% of ACDM haloes have a more prominent spatial 
-h 3D-kinematic plane than the MW one. 

In Fig. 11 we plot the properties of the rarest spatial 
-h 3D-kinematic planes that are at least as prominent as 
the MW satellite plane. In analogy to Sec. 4.5, we find that 
the planes are characterized by a large diversity in orbital 



^ACDM ( ^sat ’ ^ - ^S.s.r. ) 

Figure 12. The CDF of the incidence, /ACDM(A^sat 5 , 

^A^s.s.r.); of finding similar planes in ACDM. As in the M31 
case, this illustrates the diversity of planar configurations found 
in MW-like mocks. The MS-II and COCO results are consistent 
with the scatter expected due the low number of COCO haloes. 
The vertical dashed line and the grey area show the incidence and 
the Icr error for the MW plane of satellites. We find that 18lt|% 
of ACDM systems have an even lower frequency than the MW 
plane. 



^ACDM ( , > r|| , < ) 

Figure 13. As Fig. 12, but for the incidence, /ACDM(Afsat, ^rx, 
In this case, we compare the orbital pole dispersion, 
Agtdi between different planes and not the number of satellites 
sharing the same sense of rotation, As s r., as in Fig. 12. We find 
that 11^®% of ACDM systems have an even lower incidence than 
the MW plane. Thus, the low incidence of the MW plane is not 
a symptom of discrepancy with ACDM, as claimed by Pawlowski 
et al. (2014). 


pole dispersion, thickness and radial extent. A plane can 
be very prominent by being very thin, by having a small 
orbital pole dispersion or by a combination of the two, which 
explains the large scatter seen in the A^td and r± properties. 
Interestingly, we find that most of such planes (43%) have 
Aaat = 11, the same as the number of members in the MW 
plane of satellites whose characteristics are shown as a large 
triangle. 
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5.3 The incidence of MW-like planes 

We follow the same approach as in Sec. 4.5 and compute 
the incidence of the rarest plane around each ACDM halo. 
The outcome is presented in Fig. 12 which shows that most 
systems have planar configurations that are very infrequent, 
with over half of the haloes having a plane with an inci¬ 
dence of 0.03 or lower. The slight disagreement between the 
MS-II and the COCO results is consistent with the scatter 
expected for the much smaller sample of COCO systems. 
This appears as a systematic shift due to the correlations 
between points in the CDF. For this test, the MW plane of 
satellites has a frequency of 10“^ (vertical dashed 

line in Fig. 12) that is consistent with ACDM expectations, 
since 18l|% of haloes have an even lower incidence. 

Since for the Galactic satellites we have full kinematical 
data, we can define a new incidence, /ACDM(Afsat, ^ry, 
^Astd), that includes the 3D kinematics. This is similar to 
the incidence introduced in Sec. 4.5, except that now we 
compare Astd between different planes instead of W.s.r.- The 
CDF of this new frequency is shown in Fig. 13. Compared to 
Fig. 12, the outcome is very similar except for a slight shift 
in the CDF towards the left, i.e. towards lower frequencies. 

The new statistics indicates that the MW plane of satel¬ 
lites has an incidence of l.hlj ® ^ 10“®, as shown by the 
vertical dashed line in Fig. 13. This value agrees with the 
result of Pawlowski et al. (2014, entry 12 of their Table 5), 
whose test is very similar to ours, except that those authors 
considered the orbital pole dispersion of only 8 out of 11 
satellites. Pawlowski et al claimed that this low incidence 
of the MW plane of satellites is indicative of a shortcoming 
of the ACDM paradigm. Instead, we find that the low in¬ 
cidence indicates that planes of satellites are very diverse. 
The distribution of satellites around the MW agrees with 
ACDM predictions since llltf% of similar mass haloes have 
an even lower frequency. 


6 DISCUSSION AND CONCLUSIONS 

In this paper, we have investigated the prevalence and prop¬ 
erties of planar configurations of satellites around galac¬ 
tic mass haloes. Using two very high resolution cosmolog¬ 
ical simulations, MS-II and COCO, we have built mock 
galaxy catalogues corresponding to the satellite distribu¬ 
tions around the MW and M31. MW-like mocks consist of 
the most massive 11 satellites found outside an assumed zone 
of avoidance, mimicking the observations of the 11 classical 
Galactic satellites. M31-like mocks are modelled according 
to the PAndAS footprint and consist of the most massive 
27 satellites found within that region, corresponding to the 
M31 satellites observed by PAndAS. 

We identify the most prominent planar configuration 
as the subsample of satellites whose spatial and kinematical 
distribution is the least likely to be a statistical fluctua¬ 
tion. Applying our approach to the MW and M31 observa¬ 
tions results in the same planar distributions as determined 
by Kroupa, Theis & Boily (2005) and Ibatal3, respectively, 
even though those studies used different identification meth¬ 
ods. The good agreement is possibly due to a posteriori 
selection bias, since those authors may have inadvertently 
tuned their methods to maximize the significance of the de¬ 


tection. This would result in all the methods converging to 
the same planes. 

We have found that planar configurations of satellites 
are very common around ACDM haloes, and, moreover, ap¬ 
proximately 5 and 9% of haloes have even more prominent 
planes than those found in the MW and M31, respectively 
(see Table 1, and Figs. 3 and 10). The look-elsewhere effect 
is crucial in assessing the detection significance of a planar 
distribution, i.e. in estimating the probability of obtaining 
such a structure in an isotropic distribution (Sec. 4.3). By 
neglecting this effect, one can easily overestimate the signif¬ 
icance level by factors of ~30 and ~100 for the MW and 
M31 planes respectively. For example, while the M31 plane 
was originally reported to have a 99.998% significance (4.3a 
detection; Ibatal3), accounting for the look-elsewhere effect 
results in a more modest 99.7% significance (2.9 ct detection). 

While ubiquitous, the planes of satellites show a large 
diversity in characteristics, e.g. in the number of members, 
the plane thickness and radial extent, as well as the kine¬ 
matical structure (see Figs. 5 and 11). Most planar con¬ 
figurations are distinct, which has two major implications. 
Firstly, the notion of a representative plane of satellites does 
not exist since one cannot find a majority of ACDM haloes 
that have the same planar configuration. Secondly, the large 
diversity of planes precludes using one or two observed sys¬ 
tems for testing the cosmological paradigm on small scales. 
For such a test, a large sample of satellite systems would 
be needed to obtain a statistical measure of the system-to- 
system variation. 

The diversity of the planes of satellites is also the root 
cause behind previous claims that planes found in obser¬ 
vations are inconsistent with ACDM (Ibata et al. 2014c; 
Pawlowski et al. 2014). These authors computed the inci¬ 
dence of the MW and M31 plane of satellites to find out 
that roughly only 1 out of 1000 ACDM systems have such 
planes. To understand this result, we have computed the 
incidence of the rarest plane of satellites identified around 
each ACDM halo. We have found that the majority of pla¬ 
nar configurations have a very low incidence and that 11 and 
5% of ACDM planes have even lower incidence than that of 
the planes of satellite in the MW and M31, respectively (see 
Table 2, and Figs. 6 and 13). The low incidence is a manifes¬ 
tation of the diversity of satellite planes and, thus, contrary 
to the claim by Pawlowski et al. (2014), it does not rule out 
the ACDM paradigm. 

While the planes of satellites around the MW and M31 
are consistent with ACDM, both systems fall in the 10% tail 
of the distribution. If both planes were independent of each 
other, one might argue that there is only a ~1% chance that 
both systems are randomly drawn from a ACDM distribu¬ 
tion. This interpretation is problematic for at least two rea¬ 
sons. Firstly, both the MW and the M31 are located in the 
same large-scale environment, which in turn determines the 
preferential directions of satellite accretion (Libeskind et al. 
2014, 2015). Thus, if the environment is especially conducive 
to the formation of prominent satellite planes, then it may 
not be surprising that both systems host prominent planes. 
Secondly, the tests used to assess the prominence of these 
planes were designed a posteriori, after investigating the ob¬ 
servational data. This is especially true for the definition of 
the orbital pole dispersion that has been motivated by ex¬ 
amining the MW data. Such an a posteriori approach incurs 
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the danger of designing tests that are specifically matched 
to the peculiarities of a particular system and are not char¬ 
acteristic of the population as a whole. 

Our analysis has shown that the planes of satellites 
identihed in the MW and M31 are consistent with ACDM 
predictions based on high resolution cosmological simula¬ 
tions. This agrees with the results of Cautun et al. (2015), 
which compared the spatial and kinematical distributions of 
satellites around a large sample of isolated galaxies in SDSS 
to find agreement between observations and theoretical pre¬ 
dictions. Previous claims of an inconsistency with ACDM 
are based on misinterpreting the low incidence of satellite 
planes (e.g. Ibata et al. 2014c; Pawlowski et al. 2014) and 
on non-robust detections (Ibata et al. 2014a, see Sec. 1 and 
Cautun et al. 2015 for details). Thus, there is no convinc¬ 
ing evidence for a discrepancy between observed planes of 
satellites and the ACDM predictions. 
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APPENDIX A: PLANE IDENTIFICATION 

Here we present the practical implementation of the plane 
identification procedure. 

Al Selecting subsets of satellites 

We first describe how we identify the interesting subsets of 
satellites, which, in the next step, are used to find the rarest 
planar configurations. The simplest approach would be to 
take into acconnt every possible combination of Nsat satel¬ 
lites out of a maximum of Amax objects, with 3 Nsat ^ 
Amax- Planes with 2 or fewer objects are not physically in¬ 
teresting since any two satellites will determine a plane of 
thickness r± = 0. This naive approach, however, would re¬ 
sult in a very large number of combinations that need to be 
considered, since for fixed Asat the number of unique com¬ 
binations is 


Asat!(An,ax- Asat)! ■ ^ ’ 

In the case of the M31 system, we have Amax = 27 satellites, 
so choosing Asat = 14 would result in 2 x 10^ subsets that 
need to be considered. This analysis would have to be done 
for many thousands of systems, and for each we would need 
to generate 10® isotropic distributions. Such an approach is 
not feasible in practice. 

To overcome the immense computational challenge de¬ 
scribed above, we consider only configurations in which the 
plane members are the closest satellites to the plane. Thus, 
no other galaxy can be found closer to the plane than the 
furthest plane member. This is in line with the plane defini¬ 
tions used in earlier studies (e.g. Bahl & Baumgardt 2014; 
Gillet et al. 2015). We start by selecting a sample of A 


planes centred on the host galaxy and characterized by the 
normal vector, ripiane. To obtain these planes, we gener¬ 
ate normal vectors that are uniformly distributed on half a 
sphere, since the opposite hemisphere corresponds to iden¬ 
tical planes flipped upside down. For each such plane, we 
order the Amax satellites according to their distance to the 
plane. The interesting subsets of satellites are those made of 
the closest 3, 4, ..., Amax objects from each plane. 

To make sure that we identify all possible satellite sub¬ 
sets, we would like to have a very large number of random 
planes, A. In turn, increasing A incurs a significantly larger 
computational cost and ends up adding mostly duplicate 
subsets of satellites, which were already identified for small 
values of A. We found the best compromise to be A = 10®, 
which is large enough to contain a significant fraction of all 
possible subsets. Using A = 10® we find 93% (70%) of the 
subsets we would identify using A = 10"^ for Amax = 11 
(27), which corresponds to the total number of satellites in 
the MW (M31) system. This means that for some systems we 
are missing the satellite subset corresponding to the rarest 
plane. In those cases we end up identifying the second or the 
third rarest planes as the most prominent planar configura¬ 
tions. Using a small sample of around 200 ACDM haloes, 
we have checked that using A = 10^ instead of A = 10® 
brings only minor changes to the GDF of the prominence, 
J,rarest^ of the rarest plane (e.g. Fig. 2) and to the GDF of 
the frequency, /acdm, of those planes (e.g. Fig. 6). Thus, 
any missing subsets of satellites will not change our overall 
conclusion. 

The subsets of satellites used to compute the frequency 
of the rarest planes (Secs. 4.5 and 5.3) were identified em¬ 
ploying the same procedure except that we used A = 10®. 
That is, we used 10® random planes uniformly distributed 
on a hemisphere, which is the same as the approach used by 
Bahl & Baumgardt (2014) and subsequent studies. 

A2 Generating isotropic distributions 

Each isotropic realization is generated by picking random 
polar and azimuthal angles® for each satellite, while keeping 
constant the radial distance from the host galaxy®. Thus, 
each isotropic realization has the same radial distribution 
of satellites as the original system. This point is crucial, 
since the radial distribution of satellites has a strong effect 
on the thickness of the resulting planes. Radially concen¬ 
trated satellite distributions result in thinner planes than 
more radially extended ones. Thus, we need to generate new 
isotropic realizations for each system of satellites. 

When constructing the isotropic distributions, we also 
generate random 3D velocities, which are used for comput¬ 
ing the distribution of orbital pole dispersions, Agtd. Since 
we are only interested in the direction of the orbital momen¬ 
tum, the magnitude of the velocity is not important. Thus, 
the velocities are generated by picking random polar and 

^ The cosine of the polar angle and the azimuthal angle are se¬ 
lected from a uniform distribution spanning the interval [—1,1] 
and [0, 27r], respectively. 

® Each random realization lies within the survey mask, which is 
PAndAS for the M31 and a 19.5° obscuration angle for the MW. 
If a random point falls outside the mask, we generate new random 
angles till the point is located within the survey mask. 
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Figure Al. The PDF of the plane thickness, for isotropic 
satellite distributions inside the PAndAS survey footprint. We 
show results for planar configurations that contain Nsa,t = 8, 15 
and 27 satellites out of a maximum of Amax = 27 satellites. The 
vertical arrow and shaded region illustrate the probability that a 
fictitious plane A, which has Agat = 15 and r±_ = 17.5 kpc, is due 
to a statistical fluctuation (see the text for details). 


azimuthal angles for each satellite, with the two angles fully 
independent from the random polar and azimuthal angles 
used to obtain the position of each satellite. 


A3 The probability distribution of statistical 
fluctuations 

We now describe how to estimate the probability that the 
spatial or kinematical distribution of a set of Agat satellites 
is the result of a statistical fluctuation. This probability is 
computed using isotropic distributions, which characterize 
the degree of planarity expected from chance alignments and 
from the discreteness of the satellite distribution. 

The probability of obtaining by chance a plane of Agat 
that is thinner than r± is given by 

Agat) = r (rV) drl , (A2) 

Jo 

where the integrand is the PDF of obtaining in an isotropic 
distribution planes with Agat members and thickness, r'^. To 
compute the PDF, for each halo we generate 10® isotropic 
realizations using the procedure described in Appendix A2. 
For each such realization we find the thinnest plane with 
Agat members. The corresponding histogram over all real¬ 
izations gives the PDF of r'j^ values. The resulting PDF, for 
the case of the M31 system (Amax = 27 satellites), is shown 
in Fig. Al. For clarity, we only give the planes with Agat = 8, 
15 and 27 members. The figure also illustrates, in an intu¬ 
itive fashion, the meaning of Eq. (A2). We exemplify this 
using a fictitious plane A that contains Agat = 15 members 
and whose thickness is shown with a vertical solid arrow. 
Eq. (A2) corresponds to the shaded area to the left of the 
solid arrow. 

The probability of obtaining by chance a configuration 
of Agat satellites in which at least Na.s.i. members share the 
same sense of rotation is given by the binomial distribution 


Figure A2. The PDF of the angular dispersion of orbital poles, 
Agtd, for an isotropic satellite distribution. We show planar con¬ 
figurations that contain Agat = 5, 8 and 11 satellites out of a 
maximum of Amax = 11 satellites. To obtain these results, we 
used the same survey footprint and number of satellites as em¬ 
ployed for the analysis of the MW satellite system in Sec. 5. 


with a success probability of 0.5 . Thus, 

Ag.gm I Agat) = 2 ^ 


Ag. 


.!(Agat-Ag.g.,.)! 


(A3) 


where the first factor of 2 comes from the fact that we do not 
fix a preferential sense of rotation, allowing both clockwise 
and counter-clockwise rotations. 

The probability of obtaining by chance a plane with 
Agat members that has an orbital pole dispersion less than 
Agtd is calculated as 


r^std 

P ry I Agat) = / PDF^T-Z-^a.t (^Btd) dAltd ■ 

(A4) 

The integrand gives the PDF of the orbital pole dispersion, 
A't^j, for an isotropic distribution, which is estimated us¬ 
ing 10® random realizations, as for Eq. (A2). For each such 
isotropic realization, we find the subsample of Agat satellites 
that has the lowest orbital pole dispersion. The histogram 
of the lowest A'td values over all realizations gives the PDF 
used in Eq. (A4). In Fig. A2 we illustrate the outcome of 
such a calculation for the case of Agat = 5, 8 and 11 out of 
a maximum satellite count, Amax = 11. 
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